function [pL,qL,pR,qR] = bond_pricing_pde_boundary_conditions(w0,v_w0,w1,v_w1,lambda_input, par_bonds, reversal)
% Provide the needed boundary conditions.
if(nargin<7)
   reversal = false;   % note: if reversal = true, then we reverse the definition of lambda as lambda_star - lambda, and solve backward. 
end

qL = 0;  qR = 0;

lambda_star = par_bonds.lambda_star;
if(reversal)
    lambda = lambda_star - lambda_input;
else
    lambda = lambda_input;
end

pL =  v_w0 - par_bonds.v_w0_fun( lambda ) ;    
pR =  v_w1 - par_bonds.v_w1_fun( lambda ) ;    

% Note: the left boundary is w=0 (for numerical reasons, we have w=0.001), and the right boundary is w=1 (for numerical
% reasons, we have w=0.99);


end